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Abstract 

We describe novel CHRDU mutations in ten families with X-linked megalocornea (IVIGCI). Our mutation-positive cohort 
enabled us to establish ultrasonography as a reliable clinical diagnostic tool to distinguish between MGCl and primary 
congenital glaucoma (PCG). IVlegalocornea is also a feature of Neuhauser or megalocornea-mental retardation {MMR) 
syndrome, a rare condition of unknown etiology. In a male patient diagnosed with IVIIVIR, we performed targeted and whole 
exome sequencing (WES) and identified a novel missense mutation in CHRDU that accounts for his IVIGCI phenotype but 
not his non-ocular features. This finding suggests that MMR syndrome, in some cases, may be di- or multigenic. MGCl 
patients have reduced central corneal thickness (CCT); however no X-linked loci have been associated with CCT, possibly 
because the majority of genome-wide association studies (GWAS) overlook the X-chromosome. We therefore explored 
whether variants on the X-chromosome are associated with CCT. We found rsl 499563 16, in intron 6 of CHRDU, to be the 
most significantly associated single nucleotide polymorphism (SNP) (p = 6.81 xlO~^) on the X-chromosome. However, this 
association was not replicated in a smaller subset of whole genome sequenced samples. This study highlights the 
importance of including X-chromosome SNP data in GWAS to identify potential loci associated with quantitative traits or 
disease risk. 
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Introduction 

X-linked megalocornea (MGCl; MIM 309300) is an inherited 
congenital disorder, characterised by bilateral enlarged corneas 
with a horizontal diameter of & 1 3 mm (measured after the age of 
two years) and reduced central corneal thickness in the absence of 
raised intraocular pressure (lOP) [1] . Adult-onset cataract typically 
develops between 30-50 years [2]. Other secondary changes 
include mosaic corneal degeneration (shagreen), corneal arcus 
juvenilis, lens dislocation and mUd iris atrophy with pigment 



dispersion [3]. The condition was genetically linked to the long 
arm of the X-chromosome over twenty years ago (Xql2-q26; 
MGCl) [2,4] but the underlying genetic cause, mutations in 
CHRDU (MIM 300350), has only recendy been discovered [3]. 
CHRDLl encodes ventroptin (or neuralin-1, neurogenesin- 1 , 
chordin-like 1), a secreted bone morphogenetic protein (BMP) 
antagonist [5]. To date, presumed loss-of-function mutations in 
CHRDLl have been described in eight unrelated families affected 
with MGCl [3,6]. 
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Distinguishing megalocomea from primary congenital glauco- 
ma (PCG) in infants is clinically challenging due to overlapping 
phenotypic features. PCG is a genetically heterogeneous condi- 
tion, often associated with a broad spectrum of additional 
phenotypes [7]. PCG is characterised by high lOP and enlarged 
globe size (buphthalmos) resulting from an obstruction of aqueous 
outflow from the anterior segment. Late detection and delayed 
treatment lead to irreversible and progressive damage to the optic 
nerve and permanent visual loss can occur [8] . Genetic analysis of 
the condition has resulted in identification of three PCG loci on 
chromosomes 2p22 (GLC3A, MIM 613086) [9], lp36 (GLC3B, 
MIM 600975) [10] and 14q24.3-q31.1 (GLC3C, GLC3D, MIM 
613085, MIM 613086) [11]. Recessive mutations in CYPIBI 
(GLC3A) are the most common cause, especially in populations 
where consanguinity is common, such as the Saudi population 
[12]. Variants in MYOC have also been associated with PCG 
either independendy or in association with CYPIBI [13]. 
However, variants in this gene are more commonly associated 
with juvenile and adult-onset primary open angle glaucoma 
(GLCIA, MIM 137750) [14]. Recessive mutations in LTBP2 
(GLC3D, MIM 613086) have been identified as a cause of early- 
onset glaucoma secondary to congenital megalocornea in multiple 
consanguineous famihes [15-17]. Notably, in the eight reported 
unrelated MGC 1 families with CHRDLl mutations, none of the 
affected males (ranging from 8-72 years of age) had raised lOP, 
glaucoma or significant visual loss [3,6] . 

Megalocornea is a defining feature of megalocornea-mental 
retardation (MMR) syndrome [MIM 249310], also known as 
Neuhauser syndrome [18]. The condition is rare and characterised 
by megalocornea in addition to neurological symptoms including 
intellectual disability (ID), hypotonia and seizures [1,18,19]. 
Syndromic facial features are also described and include a 
prominent forehead, broad nasal root, epicanthus, large low set 
ears, long upper lip and anteverted nostrils [1]. The genetic cause 
of MMR syndrome is currently unknown. The majority of 
reported cases are consistent with an autosomal recessive mode of 
inheritance or a de novo cause, and the phenotypic heterogeneity 
associated with cases reported to date suggests that the condition 
may be genetically heterogeneous [19]. 

Reduced central corneal thickness (CCT) is a common feature 
of all MGCl patients with a CHRDLl mutation [3,6]. In the 
general population, CCT is known to be a highly heritable, 
normally distributed quantitative trait, and is the focus of many 
genome-wide association studies (GWAS) because it is known to be 
a major risk factor and endophenotype for open-angle glaucoma, 
one of the leading causes of irreversible blindness worldwide [20— 
25] . The X-chromosome is commonly excluded from such GWAS 
analysis [26] and given that reduced CCT is a key feature of 
MGCl, here we tested whether any common variants on the X- 
chromosome influence CCT in the general population. 

Materials and Methods 

Study subjects and clinical examination 

All investigations were conducted in accordance with the 
principles of the Declaration of Helsinki. The study was approved 
by the following local ethics committees; Moorfields and 
Whittington Hospital UK, The Royal Victorian Eye and Ear 
Hospital Australia, Odense University Hospital Denmark and 
University Hospital of Wales UK. Informed written consent, 
including permission to publish photographs, was obtained from 
all participating individuals or parental guardians on behalf of the 
minors enrolled in this study. Blood samples were donated and 
genomic DNA was extracted from peripheral blood lymphocytes 



using conventional methodologies. Patients and their relatives 
were clinically assessed by experienced ophthalmologists. Standard 
evaluation consisted of detailed ophthalmic examination and the 

additional measurement of the axial length of the eye and imaging 
of the anterior segment of the eye performed with ocular 
coherence tomography (OCT; Visante, Carl Zeiss Meditec), b- 
scan ultrasonography, and optical interferometry (lOLMaster, 
Carl Zeiss Meditec). Proband of Family K was clinically assessed 
by a geneticist (DTP), a paediatric ophthalmologist (DL) and a 
paediatrician (MG). He had standard MRI with a 1.5 Tesla 
scanner. 

Sanger sequencing 

CHRDLl was direcdy Sanger sequenced from PGR amplimers 
as previously described [3]. Segregation analyses of disease- 
associated variants in CHRDLl were tested in additional family 
members, where available, by direct sequencing of the specific 
exons carrying the mutation. The extent of the deletions 
encompassing CHRDLl was refined by targeted proximal and 
distal PCRs amplifying the CHRDLl flanking genes RGAGl and 
PAK3 (Families I and J). Exon 4 of the proximal gene RGAGl was 
ampKfied witii primers F: 5'-AAGGGTGAAGGCAACAAGG-3' 
and R: 5'-GCAAAGTGTCTTGATCTGCTAAG-3', and exon 1 
of the distal gene PAK3 was amplified with primers F: 5'- 
AGCAGAGAAGGGCTAGGGAG-3' and R: 5'-GTCTAGGG- 
\TTTGACCAAGCG-3'. To confirm the deletion identified in 
Family H, a long range PCR was performed using BIO-X-ACT 
Long DNA Polymerase (Bioline) according to the manufacturer's 
recommendations with primers F: 5'-TGAAGCCTAGAGATG- 
CAAAGTG-3' in intron 4 and R: 5'-GATGGCCACAGCT- 
CAGTCTA-3' in intron 5. To determine the breakpoint positions, 
the resulting long range PCR product was directiy Sanger 
sequenced. CHRDLl cDNA is numbered in accordance with 
Ensembl transcript ID ENST00000372042, with -1-1 correspond- 
ing to the A of the ATG translation initiation codon in the 
reference sequences. 

Whole Exome Sequencing (WES) 

WES was performed using lUumina TruSeq exome enrichment 
(v3) and a HiSeq2000 sequencer (lUumina). Reads were aligned to 
the hgl9 human reference serjuence using Novoalign (Novocraft, 
www.novocraft.com) version 2.05. The ANNOVAR tool (Open- 
Bioinformatics) was used to annotate SNPs and small insertions/ 
deletions. Filtering was performed to identify variants with a minor 
allele frequency (MAF)£0.01 in 1000 Genomes Project (www. 
1000genomes.org/), 6500 NHLBI Exome Sequencing Project 
(http://evs.gs.washington.edu/EVS/), and our internal control 
dataset consisting of 304 exomes. Variants were then prioritized 
based on their effects on the protein function (frameshift, non- 
synonymous, splice site, and synonymous). WES analysis was also 
extended to include 5'UTR and 3'UTR \ ariaiits, where txjvered. 
Variants were then cross-referenced with genes that have been 
reported to be associated with ID and/or hypotonia and/or 
seizures and/or epilepsy. The table of candidate genes was 
generated using the KEGG disease database (www.genomejp/ 
kegg/) and with reference to the report by Piton el al. [2 7] . All rare 
variants identified by WES in Individual 111:4 (Family K) in genes 
associated with th(; neurological phenotypes were verified by PCR 
amplification and Sanger sequencing using standard methodolo- 
gies (primer sequences are available on request). Segregation of the 
rare variants was tested in the proband's mother (Individual 11:2) 
by the same method (Table SI in File SI). 
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In silico analysis 

Non-synonymous missense variants identified were scored for 
liiieiy pathogenicity using SIFT (http://sift.jcvi.org/), PoIyPhen2 
(http://genetics.bwh.harvard.edu/pph2/) and Blosum62. 

Association study 

To determine whether any common variants on the X- 
chromosome influence CCT in the general population, we 
analysed SNP array genotyping data from the TwinsUK cohort 
[28]. Subjects were genotyped in two different batches of 
approximately the same size, using two genotyping platforms 
from lUumina: 300K Duo and HumanHap6 1 0-Quad arrays. 
CCT was measured in this cohort using ultrasound pachymetry 
and recorded for both eyes. Measurements were performed using 
a DGH Technology (model 500). Data collected from 1,957 
female Caucasian participants were analysed. Whole genome 
imputation of the genotypes (2,729 SNPs) was performed using 
lOOOG haplotypes with the IMPUTE v2 software. Stringent 
quality control (QC) measures were implemented prior the 
imputation, including minimum genotyping success rate (>95%) 
and Hardy- Weinberg equilibrium (P>10 % Regression analysis 
for association was performed using mach2qd and the results were 
adjusted for age effects and family relatedness among the samples. 

We checked whether this association could be replicated in a 
subset of directly whole genome sequenced individuals, made 
available through the UKIOK project (http://www.uklOk.org/). 
A total of 760 female participants with both sequencing and 
plu'not)'jMng data were analysed. A regression model for 
association was fitted using PLINK (http://pngu.mgh. harvard. 
edu/~purcell/plink/) using as predictor each single variant and as 
outcome CCT adjusted for age effects. 

Results 

MGC1 families 

We ascertained ten new families (A-J) with a diagnosis of 
MGCl (Figure 1; Table 1), the largest cohort reported to date. 
Ophthalmic features of all patients are described in Table 1. 
Consistent with our previous findings [3], all patients were 
identified as having large horizontal corneal diameters (ranging 
from 13-16 mm, white to white) and a below average CCT, 
ranging from 346-475 |im with the mean CCT being 418 |tm 
(lSD = 33nm; number eyes measured = 30) (Table 1 and [3]). 
The reported mean value for CCT in normal male and female 
adult eyes is 535 |xm (ISD = 31 |J,m) [29]. Deep anterior chamber 
depths were noted in all examined MGC 1 subjects, ranging from 
4.30-6.50 mm (Table 1). In the normal population of predomi- 
nandy European children (mean age 6.7 years), the mean axial 
length in right eyes was 22.61 (-1-/ — 0.02) mm and the mean 
anterior chamber depth was 3.34 (-1-/ — 0.01) mm [30]. In an older 
population of white adults (>65 years) mean axial length and 
anterior chamber depth were found to be 23.69 (+/— 1.16) mm 
and 3.11 (+/-0.37) mm, respectively [31]. lOPs for all MGCl 
subjects were determined to be within normal limits. Cataract 
surgery had been performed in cdl three individuals above the age 
of 40 years. Additional secondary features of the condition 
including shagreen, arcus juvenilis, and mild iris atrophy with 
pigment dispersion were also age-related (Table 1; Figure SI in 
FUe SI). 

Novel CHRDLl mutations causing MGCl 

For each affected proband (Families A-J), CHRDLl exons and 
intron-exon boundaries were screened by bi-directional Sanger 



sequencing. In total, ten CHRDLl mutations were identified in 
each of the ten respective families investigated (Figure 1; Table 2). 

Affected Individual V:4 from Family A was identified as 
harbouring a novel missense mutation in CHRDLl , c.872G>A, 
p.(Cys291Tyr). Ventroptin has three highly conserved cysteine- 
rich von Willebrand factor, type C (VWFC) domains, and Cys291 
is within the third domain (Figures lA, 2D and S2 in File SI). A 
strong X-hnked family history of megalocomea is reported in 
Family A, however further samples were not available for 
segregation analysis. Secondary changes in adults with megalo- 
comea include an appearance similar to pigment dispersion 
syndrome, but without a Krukenberg spindle on the corneal 
endothelium. Proband (V:4) was a high myope (RE —8.00, LE — 
9.00), who exhibited unusually marked posterior iris bowing with 
atrophy, peripheral iris transillumination defects and increased 
trabecular meshwork pigmentation (Figure SI in File SI). Aged 43 
years he had bilateral phacodonesis with nuclear sclerotic 
cataracts, but normal pressures and no evidence of glaucomatous 
optic neuropathy. 

The proband in Family B (Individual IV:5) was found to have a 
unique hemizygous nonsense mutation in CHRDLl, c.240T>A, 
p.(Cys80*). Segregation analysis revealed that the proband's 
mother (111:7), maternal grandmother (11:5) and maternal half 
great aunt (111:2) were all heterozygous carriers of this mutation. 
His sister (IV: 4) was found to be wild- type for the variant. 
Interestingly, affected relatives of this Danish family were 
originally reported in the 1960s [32]. The proband's affected 
great maternal grandfather, Individual 1:2 (Figure lA), corre- 
sponds with Individual lV-7 in the original article with a greatiy 
extended version of the pedigree spanning 7 generations and 
detailing many affected males [32]. 

In Family C, Individual IILI and his maternal grandfather (1:1) 
were found to have a novel 2 bp deletion affecting the splice 
acceptor site of CHRDLl exon 12; c.l247-l_1247del. This 
mutation is predicted to cause aberrant pre-mRNA splicing of the 
CHRDLl transcript in vivo. His two younger brothers have also 
been diagnosed with the same condition; however their DNA 
samples were not available for testing (Figure lA). The probands 
in Family D (Individual 11:1) and Family E (Individual 11:1) were 
both identified to have novel hemizygous nonsense mutations in 
CHRDLl; c.l00G>T, p.(Glu34*) and C.2970A, p.(Cys99*), 
respectively (Figure lA). Individual 11:1 in Family F was found to 
have a unique missense mutation in CHRDLl, c.865T>C, 
p.(Cys289Arg) (Figure lA). Similar to other missense mutations 
(Family A and [3]), the affected residue is positioned within a 
highly conserved cysteine-rich V'WFC (Figure 2 and S2 in File SI). 
A novel nonsense mutation in CHRDLl, C.2290T, p.(Arg77*) 
was identified in Individual 11:2 in Family G. The proband's 
maternal uncle was reported to be affected with megalocomea, 
however, no DNA was available for testing. 

The proband in Family H (Individual 11:1) was initially 
identified as harbouring a deletion of (CHRDLl exon 5 by PCR. 
In order to confirm the deletion and define the deletion 
breakpoints, a fragment of DNA spanning intron 4 to intron 5 
was amplified and sequenced, revealing a 9,033 bp deletion and 
8 bp insertion; g.l09964279_109973310delinsTTATGCTA. The 
deletion encompassed exon 5, and partial sequence from intron 4 
and intron 5, and is predicted to result in a frameshift and the 
introduction of premature termination codon, p.(Glul01Glyfs*42) 
(Figure IB). Segregation analysis of the deletion confirmed that the 
proband's mother is a carrier and the father is wild-type. 

Larger deletions encompassing the entire CHRDLl gene (12 
exons) were identified in the affected males in Families I and J, by 
PCR. Successful amplification of the closest exons of neighbouring 
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Figure 1. X-linked Megalocornea Families A-J. (A) Pedigree of Families A-G. Black shaded squares denote clinically and genetically confirmed 
affected males; grey shaded squares denote clinically diagnosed affected males but DNA samples were not available for testing; dotted circles denote 
genetically confirmed carrier females; ? = presumed carrier females but DNA samples were not available for testing; clear squares and circles denote 
unaffected individuals. Arrowhead indicates proband in each family. Control sequence electropherogram is shown above patient sequence. (B) 
Pedigree of Family H and sequence electropherogram showing a 9,033 bp deletion encompassing CHRDLl exon 5 and 8 bp insertion. (C) Pedigree of 
Families I and J. Deletion of the entire CHRDLl gene (exon 1 to exon 12) in the proband is shown in agarose gel images. The flanking genes, RGAGl 
and PAK3 are present. NTC, non-template control (D) Schematic representation of presence or absence of the CHRDLl gene and flanking genes 
RGAGl and PAK3 (EnsembI nomenclature hg19 genome build) in Families I and J. 
doi:10.1371/journal.pone.0104163.g001 
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Figure 2. Summary of CHRDL 1 mutations. (A) Schematic of CHRDLI protein domains. The following abbreviations are used: SP, signal peptide; 
VWFC, von Willebrand factor, type C domain. (B) Schematic of the CHRDLI gene. (C) CHRDLI mutations previously reported in X-linked megalocornea 
families [3,6]. Frameshift, splicing, nonsense, missense, and whole gene deletion mutations were identified. (D) Novel CHRDLI mutations identified in 
MGC1 families in this study (Families A-J). The dagger (t) indicates Family K with iVlMR syndrome. 
doi:1 0.1 371/journal.pone.01 041 63.g002 



genes, RGAGl and PAK3, demonstrated in both instances that 
the deletions identified do not encompass any additional currently 
annotated genes (Figure IC and D). 

Ultrasonography as a diagnostic tool for MGC1 

Our molecular diagnosis of MGC 1 patients, here and previously 
[3] , has allowed us to define the key phenotypic characteristics of 
MGCl (Table 1) and to difierentiate them from PCG (Figxire 3). 
Patients with MGCl have a large corneal diameter and thin 
cornea, but without corneal oedema or evidence of breaks in 
Descemet's layer. The depth of the anterior chamber in MGC 1 is 
also typically significantly greater than in patients with PCG. 
Ultrasonography is a reliable clinical diagnostic tool to distinguish 
the two conditions and, if on ultrasonography, the ratio of the 
anterior chamber depth to the total axial length is SO. 19 then a 
diagnosis of MGCl is extremely likely, unless there is gross 
coexisting axial myopia (>36 mm) (Figure 3; Table 1). 

MMR syndrome 

The male proband in Family K is the only child of his parents, 
who are unrelated. Pregnancy was compKcated by 7 days of 
vaginal bleeding at 8 weeks gestation, and persistent abdominal 
pain from 1 7 weeks of pregnancy. Due to the abdominal pain, 
labour was induced at 39 weeks gestation. Delivery was normal 
and he weighed 3.917 kg (91'"' centile). In the first week of life, he 
had some choking episodes, which were thought to be due to 
gastro-oesophageal reflux. Thereafter, feeding and weight gain 
were satisfactory. He had an apnoeic attack at the age of 7 weeks. 
At the age of 3 months, he developed epilepsy and was started on 



Sodium Valproate. By 6 months of age, there were concerns about 
his psychomotor development. He was found to have global 
developmental delay; he only walked at the age of 20 months and 
he required speech therapy. Bilateral megalocorneae were 
diagnosed in infancy. Physical features at the age of 3 years 
included a broad forehead, bilateral megalocorneae, bilateral 
epicanthic folds, a tented upper lip and downturned corners of the 
mouth (Figure 4B). His occipital frontal circumference (OFC) was 
on the 25''^ centile. Magnetic resonance imaging (MRI) of his 
brain and karyotying were reported as normal. He was also very 
hypotonic initially, but this has improved. On the basis of his 
medical history and physical features, a diagnosis of MMR 
syndrome was made. At the age of 10 years, he has moderate 
intellectual disability (ID) and a diagnosis of autistic spectrum 
disorder. His antiepUeptic medication had been discontinued at 
the age of 9 years, and he has had no further seizures. He has 
episodes of challenging behaviour. Growth parameters were 25*— 
50* centile for height, 9r'-98* centile for weight, and his OFC 
was between the 9* and 25* centile. 

Neither of the proband's parents have a history of megalocor- 
nea, ID or epilepsy. There is no paternal contact. His mother had 
laser iridotomies at the age of 47 years for bilateral, narrow 
anterior chamber drainage angles. He has 2 maternal half-sisters 
and a maternal half-brother, who are healthy. His maternal 
grandparents are both deceased; his maternal grandfather is 
reported to have had "large eyes", cataracts and glaucoma. He 
also has 3 paternal half-siblings. Two of them are said to have 
epilepsy and ID, although this may be due to infantile meningitis 
in one of the children. 
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CCT: 0.47 mm 





Figure 3. Comparison between biometry of X-linl<ed megalocornea and congenital glaucoma. (A) AC, anterior chamber depth (bacl< of 
cornea to front of lens); AL, axial length (anterior corneal surface to vitreo-retinal interface at the back of the macular). Transverse B-Mode ultrasound 
scans demonstrate higher ratio of the AC to AL in megalocornea, 0.204 (left) than congenital glaucoma, 0.135 (right). (B) Ocular coherence 
tomography (OCT) images comparing megalocornea (left) and arrested congenital glaucoma (right). Although there is enlargement of the anterior 
segment in both conditions, patients with megalocornea tend to have a greater anterior chamber (AC) depth, lower central corneal thickness (CCT), 
and deeper anterior chamber angles (ACA) than congenital glaucoma patients. 
doi:10.1371/journal.pone.0104163.g003 



Identification of a CHRDLl missense mutation in a male 
patient diagnosed with MIVIR syndrome 

Given the rarity of MMR syndrome and tlie phenotypic 
variability reported in the literature [19], we hypothesised that the 
condition may not be genetically homogeneous. Since megalocor- 
nea is a primary feature of MMR syndrome, we considered that in 
this affected male (Individual 111:4 Family K, Figure 4), the 
megalocornea phenotype may be due to a mutation in CHRDLl . 
Sanger sequencing of CHRDLl revealed a unique hemizygous 
missense mutation, c.464G>A, p.(Cysl55Tyr) (Table 2) in the 
proband within a conserved cysteine-rich VWFC domain. His 
unaffected mother was identified to be a heterozygous carrier 
(Figure 4A). Despite mutations in CHRDLl being previously 
associated with abnormal brain morphology [3], it was unlikely 
that his extraocular features were caused by the CHRDLl 
mutation because no other patients with a mutation in CHRDLl 
present with developmental delay or any other of the extraocular 
phenotypes associated with MMR syndrome [3]. 



WES analysis of genes previously associated with ID, 
epilepsy, seizures and hypotonia 

Next, we considered the possibility that the extraocular 
phenotypes may be due to a mutation in an X-linked intellectual 
disability (XLID) gene in linkage disequilibrium (LD) with the 
CHRDLl mutation, or an autosomal ID gene. We therefore 
performed WES with the aim of determining if any additional 
genetic variants could be contributing to his extraocular features. 

Overall 21,210 exonic and 61,579 total genomic sequence 
alterations were identified by WES, with respect to the reference 
sequence. As causal variants are hypothesised to be rare, we 
filtered the dataset for rare variants with a MAF^O.Ol in 1000 
Genomes Project, 6500 NHLBI Exome Sequencing Project, and 
our internal control data set consisting of 304 exomes. WES 
further confirmed the CHRDLl mutation. We created a database 
of genes tiiat have been associated with ID and/ or hypotonia and/ 
or seizures and/or epilepsy, and cross-referenced the filtered rare 
variants dataset against this database. Using this approach, we 
identified ten rare variants of unknown significance in genes 
associated with these neurological phenotypes including four 
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Figure 4. Novel CHRDLl missense mutation in a patient diagnosed with MiVIR. (A) Pedigree of Family K with megalocornea-mental 
retardation (MMR) syndrome. Squares, males; circles, females; diamonds, unknown gender; shaded, affected; dotted, carrier; clear, unaffected. 
Arrowhead indicates proband. Sequence electropherograms show the CHRDLl missense mutation c.464G>A; p.(Cys155Tyr), which segregates in 
Family K. (B-C) The proband at ages 3, and 6 years, respectively, presented with bilateral megalocorneae, broad forehead, bilateral epicanthic folds, a 
tented upper lip, and downturned corners of the mouth. (D-E) Frontal and sides of the proband at age of 10 years. 
doi:1 0.1 371/journal.pone.01 041 63.g004 



heterozygous missense variants, two heterozygous synonymous 
variants and four variants present witliin non-protein coding 
regions (5' and 3' UTR regions) (Table SI in File SI). Each of the 
ten variants was validated by direct Sanger sequencing in the 
proband and his unaffected mother; paternal DNA was not 
available for testing. Two rare autosomal variants in SEMA3E 
(MIM 608166) and DDBl (MIM 600045) were identified in the 
proband and his unaffected mother in the heterozygous state, 
suggesting that they are unlikely to be causative of the extraocular 
phenotypes observed (Table SI in File SI). Five of the remaining 
eight heterozygous variants were identified in genes associated 
with autosomal recessive forms of disease {ERCC8, UROCl, 
ADCK3, ERCC6, and ERCC3) suggesting that, in isolation, each 
of these rare variants would also be unlikely to be causative. 

A unique heterozygous variant (g. 161326951 A>C) was 
identified in the 3'UTR oi GABRAl (MIM 137160, Table SI 
in File SI). Non-synonymous heterozygous variants in GABRAl 
have been associated with juvenile myoclonic epilepsy [33]. 
However, considering the rarity of critical regulatory elements in 
3'UTRs, a mutation in this region is usually tolerated. It is 
therefore unlikely that this variant would be causative. The 
proband's mother (Individual 11:2) was wild-type for this variant, 
so either the proband inherited the variant from his father 
(Individual 11:3), or it is a de novo change. The proband's father 
has no reported history of epilepsy. 



A rare heterozygous missense variant, p.(Thrl016Asn), absent 
in the proband's mother, was identified in DIP2B (MIM 6 11 379, 
Table SI in File SI). An expansion of a CGG repeat in the 5'UTR 
of DIP2B has been reported in two patients with ID (FRA12A) 
[34], however no missense variants have been associated with 
disease. 

A unique hemizygous 1 bp deletion in the 5'UTR region of 
CUL4B (MIM 300304) was identified in the proband (Table SI in 
File SI). Interestingly, a 2 Kb non-coding deletion that disrupts 
the 5'UTR of CUL4B has previously been identified to co- 
segregate with ID in a single family, resulting in loss of CUL4B 
expression, suggestive of a pathogenic effect [35]. Mutations in this 
gene in several families with XLID were initially reported in 2007 
[36,37]. InterfamUial and intrafamUial phenotypic variability was 
noted with regard to several features, including the degree of ID 
and the growth parameters (height and OFC). Additional 
publications further confirmed the association with moderate to 
severe ID and marked speech delay in all patients. Challenging 
behaviour, seizures, and a prominent lower lip were also common, 
as seen in our patient [38,39] . 

Notably, CUL4B is situated approximately 900 Kb downstream 
of CHRDLl . Since the proband's mother is also a carrier of this 
variant, it is likely that the rare variants identified in CUL4B and 
CHRDLl are in LD. This suggests that MMR in this patient may 
be a digenic X-linked trait. However, the proband's deceased 
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maternal grandfather reportedly had 'large eyes' in the absence of 
extraocular features. This implies that he may have had MGCI 
and the CUL4B variant could have occurred as a de novo event in 
the proband's mother. However, the extended family pedigree 
indicates a possible paternal transmission of ID and epilepsy in the 
proband and his affected half-siblings, which would exclude an X- 
hnked gene. 

Unfortunately, we were not able to test these hypotheses due to 
lack of familial DNA samples, so the potential pathogenicity of the 
CUL4B and other variants identified, and therefore the genetic 
cause of extraocular phenotypes in the proband, remains 
unresolved. 

Association study 

CCT is a normally distributed (Figure S3 in File SI), highly- 
heritable quantitative trait in the human population [40]. Despite 
many GWAS investigating CCT-associated loci, the X-chromo- 
some has been largely ignored. We therefore wanted to determine 
whether common variants on the X-chromosome in the general 
population are associated with CCT. To test this hypothesis we 
examined SNP array genot}ping data for 1,957 female partici- 
pants of British descent from the TwinsUK cohort. Their mean 
CCT was 545.24 jxm (±40.03 \im). The most significant CCT- 
associated SNP on the entire X-chromosome was rs 1499563 16 
(p = 6.81 X 10"'*) (Figure S4 in Ffle SI). This SNP is located within 
intron 6 of CHRDLl , and the positive effect allele (A) has an MAF 
of 0.017 in the TwinsUK dataset (0.008 in the 1000 Genomes 
dataset). To assess potential functionality of the rs 1499563 16 
variant, we compared the effect of the reference and variant allele 
using pre-mRNA splicing prediction programs including Human 
Splicing Finder (HSF) [41], NNSPLICE [42], NetGene2 [43] and 
SpliceAid [44] , but there were no predicted motif alterations or 
alternative splicing activity. 

Given the known limitations of genomic imputation for variants 
with low MAF, we attempted a rephcation of these results in a 
smaller subset of 760 individuals from the TwinsUK cohort, for 
which both CCT and genome sequence data were available. We 
did not observe a significant effect in this smaller sample and were 
unable to replicate the results from the imputation. The current 
conflicting data makes it impossible to determine if this non- 
validation is the result of insufficient power from the limited 
number of samples available, accumulation of potentially influen- 
tial imputation mistakes, or if the initial imputed observation was 
purely a type one error. 

This data initially suggested that the CHRDLl locus may be a 
quantitative trait locus for CCT, however, given the uncertainties 
during the replication in direcdy acquired genotypes (as opposed 
to imputed ones), further validation is warranted before the 
CHRDLl locus can be substantiated as influencing CCT in the 
general population. 

Discussion 

All MGC 1 families ascertained in this study were found to have 
unique CHRDLl mutations including nonsense mutations 
p.(Glu34*), p.(Arg77*), p.(Cys80*) and p.(Cys99*), missense 
mutations p.(Cys289Arg) and p.(Cys291Tyr), a frameshift muta- 
tion p.(Glul01Glyfs*42), a splice site mutation (c.l247-l_1247del) 
and deletions encompassing the entire CHRDLl transcript. To 
date, all patients recruited to our study with a diagnosis of MGC 1 
were found to have a mutation in CHRDLl [3]. 

Members of the transforming growth factor beta (TGF-P) 
superfamUy, including BMPs, are known to play important roles in 
ocular development [45,46]. The concentration of active BMP 



hgands is tightiy regulated by a collection of secreted proteins, 
including the chordin family, which act as BMP antagonists. The 
chordin proteins interact with BMPs via the conserved cysteine- 
rich VWFC domains [47]. CHRDLl encodes ventroptin, which 
contains three VWFC domains. Consistent with our previous 
findings [3], all missense mutations identified to date are located at 
highly conserved cysteine residues within these domains, which are 
predicted to abrogate interaction with BMPs. 

There is an important clinical need to differentially diagnose 
MGC 1 from PCG. In our molecularly diagnosed cohort of MGC 1 
patients, we have been able to use ocular imaging define the key 
chnical characteristics that assist in differentiating MGCI from 
PCG (Figure 3). Since we have shown that CHRDLl mutations 
do not lead to increased TOP or glaucoma (Table 1) [3], the 
repeated examinations that are required to monitor lOP (under 
anaesthesia) for PCG patients are not required for MGCI patients. 
A molecular genetic test for CHRDLl can confirm a diagnosis of 
MGCI and has important implications for clinical care. 

As described, megalocornea is also a key pathognomonic feature 
of MMR syndrome. Intriguingly, we identified another novel 
missense mutation in CHRDLl , c.464G>A; p.(Cysl55Tyr) in a 
male patient with a diagnosis consistent with MMR syndrome. 
This CHRDLl mutation accounts for the ocular phenotype, and 
although the genetic cause of the extraocular phenotypes in this 
individual currentiy remains unexplained, this is the first report 
revealing a causative mutation in an MMR case. Recruitment of 
more MMR patients and a more comprehensive delineation of 
MMR syndrome are essential in order to investigate the 
underlying genetic cause(s) of this poorly defined syndrome. Our 
findings demonstrate that CHRDLl should be considered and 
screened in males diagnosed with MMR syndrome. Further WES 
and copy number variant analysis of male and female subjects will 
help to define the contribution of X-linked and autosomal genes to 
this syndrome. 

We observe that an X-linked locus may be associated with 
CCT. In addition to MGCI, corneal thinning is associated with 
rare connective tissue disorders such as Brittie cornea syndrome 
(BCS; MIM 229200) [48] and Ehlers-Danlos syndrome (EDS; 
MIM 130000) [49]. Notably, common variants in genes associated 
with certain forms of these rare conditions, ZNF469 (BCS; MIM 
612078) and C0L5A1 (EDS; MIM 120215), have been shown by 
numerous GWAS to be associated with CCT in the general 
population [20-25]. The potential association of a CHRDLl SNP 
(rs 1499563 16) with CCT mirrors these findings, and may have not 
been detected previously because non-autosomal data have largely 
been excluded from the majority of CCT GWAS studies [26]. 
Currently, the majority of GWAS focus upon analysis of 
autosomal variants only, as analysis of sex chromosomes for 
GWAS remains technically challenging [26,50]. However, given 
the lack of replii:ation in this study, it will be important to test for 
association of X-linked lo[:i, and the CHRDLl locus, with CCT in 
male cohorts, because it is possible that phenotypic effects of X- 
Unked loci (including the CHRDLl locus) may be masked in 
females, who have two copies of the X-chromosome. In males, in 
the hemizygous state, we h)'pothesise that a greater effect on 
corneal phenotypes may be observed. Our study highlights the 
potential importance of analysing X-chromosome SNP data in 
GWAS to identify loci associated with quantitative traits or disease 
risk. 

Supporting information 

File SI File includes Figures S1-S4 and Table SI. 

(DOC) 
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